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Angular distributions in the final state of ty t] photoproduction on nucleons are considered. As a 
formal base the familiar isobar model is used in which the ir°r]N state is a product of the resonance 
decay into 7/A(1232) and 7rSn(1535) channels. One of the principal assumptions used is that in the 
actual energy region the reaction is dominated by a single resonance state. The developed formalism 
can serve as a tool for testing spin and parity of that resonance. 
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I. INTRODUCTION 



Recent measurements [T| of the total cross section for n°ri photoproduction off the proton have shown that its 
value rapidly rises with increasing energy and in the region E~ > 1.2 GeV it exceeds the cross section for single rj 
photoproduction. At the same time, already a simple analysis shows that the dynamics of this reaction is by no means 
trivial. Indeed, calculations made with the Born diagrams, Fig. [lja-f), give only about 10% of the measured cross 
section value thus pointing to the crucial role of resonances in this reaction (diagrams (g) and (h) in Fig. [I]). 

The authors of Ref. have proposed as a dominant mechanism excitation of the resonance 1)33(1 700) with a 
1— 1 1 subsequent decay into the ?yA(1232) channel. It is quite clear that if the role of the Born terms is negligible, then 
43 in order to account for the rapid rise of the cross section near threshold, one has no choice but to take a resonance 
1 1 decaying into s-wave 77 A state for which the 1)33(1700) is a very good candidate. 

However, according to the experimental results of Ref. [l], the cross section does not demonstrate a pronounced 
resonance like energy dependence. It reaches about 4 ^b at E 1 = 1.3 GeV and then does not essentially change at 
least up to Ej = 2.2 GeV. Such a smooth behavior might be due to the fact that the mass of the excited resonance 
is close or even below the threshold energy, so that the corresponding peak is folded with the increasing reaction 
I ] phase space. This situation might lead to strong uncertainties in determining the resonance parameters when fitting 
■ the cross section. Furthermore we can not exclude the case with two or more strongly overlapping resonance states, 
I/"") j having different spin and parity and thus contributing incoherently. To clarify the situation a detailed analysis of 
IO 1 partial waves becomes of special importance. Some work in this direction has already been done by the Bonn-Gatchina 
group [l| . Their method based on the formalism of Ref. [3] has allowed them to analyse the role of different partial 
waves. According to their results in the region 1.07 GeV < E 1 < 1.45 GeV the dominant contribution comes from 
CO 1 the resonances L» 33 (1700) and D 33 (1940). 

The present study deals mainly with the angular dependence of the cross section which is proposed as a sensitive 
criterion to identify the underlying production mechanism. Obviously, in the event that only one partial wave 
dominates, the corresponding angular distribution is governed by the lowest powers of sine or cosine of spherical 
angles, whereas the Born sector, where large amount of partial waves are involved, provide more or less monotonic 
angular dependence. Therefore, if we are able to separate the shape of the differential cross section associated with 
pure harmonics we can expect that not only the resonance part of the amplitude can be isolated from the background 
but that it is also possible to identify the quantum numbers of individual resonance states. 

For simplicity as a first step we consider an ideal case when in the actual energy region of energy the amplitude is 
dominated by only one resonance R and contributions of other waves can be neglected. We put the resonance mass 
to Mr = 1.8 GeV and take the total width Tr = 300 MeV. The quantum numbers, spin and parity J v , are treated as 
model parameters. Our aim is to investigate the dependence of the angular distributions on the choice of R(J V ) and 
in this way to study the signature of individual partial waves. Perhaps, the next natural move would be to identify 
R(J W ) with more or less well-established resonance states known, e.g., from the PDG listing [4|. However, within the 
quite rudimentary state of the database and the theoretical descriptions, such an approach cannot lead to a unique 
solution. 

The paper is organized as follows. First, in Sect. II we briefly describe the formalism of the sequential decay of 
the resonance R according to the scheme: R — ► baryon(| + )+ meson(0~), baryon(| + ) — > baryon(i + ) + meson(0~) 
and derive expressions for the angular distributions. Then in Sect. lIIll using the assumption of the isobar model 
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FIG. 1: Diagrams for — » -K°rjN used in the calculations. The notations A and iV* stand for P33(1232) and S'ii(1535), 
respectively. 



we calculate the angular and energy dependence within different hypotheses about spin-parity of the resonance R. 
Finally, in Sect. IV we summarize our main qualitative results. 

II. ANGULAR DISTRIBUTIONS. FORMALISM. 

We consider the process 

7 (fe) + Ni(-k) - 7T°(&) + V (q n ) + N f (p f ) , (1) 

where the 3-momenta of the particles in the overall cm. frame are given in parenthesis. In general, the resonance 
mechanism of the reaction ([1]) can be realized according to the following two schemes (throughout the paper the 
resonance S , n(1535) is denoted by N*) 

(a) : 7 iV -> R{,r) -> f]A + -> ^-qN , (2) 

(b) : 7 iV -► R(J*) -► 7T°iV* -► 7r°7?iV 

(see diagrams (g) and (h) in Fig.[T]), the relative amount of which depends on the details of the reaction dynamics. 
The first scheme with R = Z?33(1700) was considered in [2| as a main driving mechanism of the reaction (JTJ. The 
second sequence appears in Ref. [2| due to strong rjN interaction via N* excitation. 

The resonance states R considered in this paper are listed in Table U together with orbital momenta associated 
with their decay into different channels. Throughout this section we assume that the ir°r] production always proceeds 
according to the scheme (a). The scheme (b) will be included in the next section where we present our results obtained 
within the isobar model. 

When writing the resonances in the form L2T2J in Table U we took into account that only the states with isospin 
T = 3/2 decay into the 77A channel. As already noted, we do not try to identify the states R^J 7 ') with the baryon 
spectrum known from PDG [ij . But if only the quantum numbers are taken into consideration the states collected in 
TableUmay be identified with, e.g., 531(1900), P 31 (1910), £> 3 3(1700), F 3 3(1920), D 3 5(1930), and F 35 (1905). These 
resonances, except for Z?33(1700), belong to the third group and can influence the low-energy region only through 
their large widths. It is also worth noting that all the mentioned states are characterized by quite a weak nN mode 
(generally less than 20%) and therefore can intensively decay into the two-meson channels. 

The kinematics of the reaction ((TJ is presented in Fig.O We select the z-axes along the photon momentum k. The 
production plane is spanned by the momenta k and <f j; , so that (f> n = 0. The decay plane of the TtN pair is fixed by the 
momenta and pf. We denote by £1 — (6, </>) the solid pion angle in the irN rest frame. To describe the A — * irN 
decay, two types of the coordinate systems Ox' y' z' are used. In the first one, the canonical frame (also referred to as 
Adair frame, Fig.^a)), all three axes are codirectional to those of the OXYZ system, so that the system O'x'y'z' is 
deduced by the Lorentz boost determined by the vector q v +pj. In the helicity frame (Fig.^b)) the z' axis is aligned 
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TABLE I: Angular momenta associated with a decay of the resonance R(J*) into hadronic channels. In the rjA states only the 
lower of two possible values of L is assumed. The resonance S l n(1535) is denoted as N*. 



L(nN) 



L( V A) 



L(nN* 



| + (^3) 



(D 3 



along the vector q v + p and the x' axis is in the production plane. To shorten the notations we denote by the polar 
angle of the vector p\ — q n + p ' = —q v in the overall cm. system, and by Q — (9, <f>) the pion momentum in the ttN 
rest frame. The coordinate systems above will further be referred as K- and if-system respectively. Clearly, they are 
connected to each other by a rotation with the angle Q around the y' axis. 

The amplitude associated with the diagram (g) in Fig. [T] has the following form in the canonical and the helicity 
frame 

£ /X (e, n) = -±= aRA * E ^4lt^k §m a Y irrM d L ML0 (e) , (3) 

V R(JL) m 



V R(.JL) rn 

Here, C^^i - ama are the usual Clebsch-Gordan coefficients for the coupling j\ + j 2 = J. The index L is the r) + A 
orbital momentum and d^^^d'l ,) are the rotation matrices. In the second equation /ka stands for the A helicity. The 
spherical function Yi TO (f2) specifies the angular dependence of the decay A — > 7riV in the 7r7V rest frame. The index A is 
the initial state helicity. Since the momenta of the initial particles are along the Z-axes, the photon helicity is combined 
with the nucleon helicity to give the z projection of the total angular momentum J equal to A € {±1/2, ±3/2}. The 
corresponding amplitude of the transition 7^ — > R( J 71 ") is denoted by AP. Parity conservation requires that 

A« A = (-1)4+J-U«. (5) 
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The parameters a characterizing the individual resonances R contain constants and energy dependent functions 
(resonance propagators, barrier penetration factors, coupling constants, etc.), which detailed structure is irrelevant 
for further discussions. We note that the nucleon magnetic quantum numbers rrif in ((3| and (j4|) are projections of 
the final nucleon spin on different z'-axes, according to their definitions in the K- and fZ-system (see Fig. [2]). 

Using Eq. (|3|) or ((4]) together with ([5]) it is easy to verify that parity conservation leads to the following symmetry 
property 

t_i£-A(e;M) = (-ir'- A Oe ; 0,-0). (6) 

The square of the matrix element can easily be written down from ([3]) and ((4]). Since in this section our main object 
is the form of angular distributions and not the absolute value of the cross section, it is convenient to introduce the 
distribution functions, normalized to unity. In the isobar model discussed in the next section we assume that the 
only nonvanishing contribution to the amplitude comes from the sole state R(J W ), whereas other resonances can be 
neglected. Deriving the formulas below we will always adhere this somewhat oversimplified picture. 

First, we consider the distribution over the angle O of the irN system in the overall cm. frame. For this purpose 
we define the distribution function 

w(e) = ^[j2 |* m/ A(e;r>)| 2 dn, (7) 

m/A 



N 



A R 



A R 



W{&) sin9 de = 1 



(8) 



Using ([3]) for a given resonance i?(J 7r ) we obtain 

2L + 1 



W{Q) 



4(1 + a) 



' M L A=±i 



JX 

L Ml 



\X-M L f+ a 



A=±; 



(C JX , 

\ LMl § A— Mi 



>')■ 



(9) 



where the parameter a is defined as 



As is evident from Eq. (JHJ), the distribution over cosO is flat for J 77 — l/2 ± since in this case a = and 

E (c 



(10) 



LM L f A — Mi, ' 



2/(2L + 1) does not depend on Ml- The same is true for L = which we have in the state 



A=±l/2 

with r = 3/2- (see Table U). Thus 



W(Q) 



if r 



2 



(L = 0) 



(ii) 



Otherwise the shape of W(Q) is described by a polynomial of cos 9 of the order 2L. Table UT1 lists the function W(Q) 
for all six transitions considered here. As is seen from the Eq. ©, the exact form of W(Q) depends on the ratio 
a = (A^, 2 /Ay 2 ) 2 , This fact brings a model dependence into our analysis, especially if the electromagnetic amplitudes 
are poorly known. On the other hand, for J n = 5/2+ all coefficients in the expansion are positively defined 
(a > 0). As a result, the corresponding distribution will always reach its maximum at | cosO | = 1. For J n = 3/2+ 
and J v = 5/2 _ the shape of the distribution might be quite sensitive to the value of a. In such a situation indirect 
information about the spin of R can be obtained from the complexity of the angular distribution, which as mentioned 
before is fixed to 2L for J > |. 

Now we turn to the £1 dependence at fixed O. The corresponding distribution function reads 

w{@]n) = NW) EK»/A(e;^)| 2 , 

m f X 

where the factor N(Q) is determined by the normalization condition 



J w{e-,ri)dn = 1. 



(12) 



(13) 
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TABLE II: Distribution W(Q) © over the angle 9 of the ttN system in the overall cm. frame and the 9 distribution 
f(9) — 2-7T 14^(6 = 0; 9, <f>) (|22[l of pions produced by the A decay in coincidence with r] mesons at 9,, = n. 



r(L 2 



W{&) 



\ ( 5 3l) 



1 + 
2 



2 



(Ps: 



3 -(F„) 



2 

5 + 1 



1 



10(l+o) 



(7 + 3a -6(1 -a) cos 2 9) 



28(1 3 +Q) (2 + 7a + 5(4 - 5a) cos 2 9 - 10(2 - 3a) cos 4 9) 
20^(2 + 3a + (4 + a) cos 2 9) 



±(l + 3cos 2 6») 
i(l + 3cos 2 6») 
3^(1 + 3a + 3(1 -a) cos 2 9) 
4(1^(1 + 27a + 3(1 - 9a) cos 2 9) 
4^1^(1 + 18a + 3(1 - 6a) cos 2 9) 
4T 3|^ y (l + 2a + (3-2a) cos 2 9) 



It is convenient to present the function W(Q; f2) in the form 

W(e;Q) = ]T Y{ m (n)p mm ,(Q)Y lm ,(Q) , (14) 

mm,' 

where, e.g., from Eq. ([3]) the correlation coefficients p mm '(0) read 

Pw(0) = Jr¥ME E M W<o(6) ^, (6) (15) 

^ ' m f XM L M' L 

with 

„M L M' L _ r 5«4 ^I-Ma r JA ^JA I R /.fll 2 

°m/A,mra' ~ 2^ °lm^m / U lm' ±m/ £M £ §Af A £M^§A^ I A I ■ U°J 

m/A 



The unit-trace condition for the matrix p mm i immediately follows from the Eqs. (| 1 3j> and (|14p 

^p rom (6) = l. (17) 

m 

It is intuitively clear that the structure of the pion angular distribution will be governed by the A spin Ja = 3/2 and 
should contain polynomials of cos6> up to the order 2Ja — 1 = 2. Using Eq. (fl"4|) one has 

W(®;0,4>) = ^- fpoocos 2 9 + -{pu + P-i-i) sin 2 9 - sin 2 (TZep^x cos 20 - Tmpx-x sin 20) 
47r V 2 

+ —= sin 26 ) 7?.e( j o_io — pio) cos0 = sin 29Tm (p-xo + pxo) sin V (18) 

V 2 V 2 ' 

where we have dropped the argument in /) mm »(6). In (|18p the hermiticity of the matrix p mm ' was already used. 
Furthermore from ^ it is evident that all elements p m m' are real and p_ TO _ m / = (_i) m -™ p mm> . Taking also into 
account the normalization condition (JX7J) we arrive at the result that from nine real elements p mm ' only three, pocb 
pio and p\-\ remain independent so that the distribution function is reduced to 

W(<d;9,<j)) = -f- fp ocos 2 6»+i(l-poo)sin 2 6'-\/2^epio sin 29 cos 20 - pi _ x sin 2 cos 2</)) . (19) 
47r V 2 / 

Projection of (JTUJ) on the y' axis gives 

^(6; 9 = = |) = A (i(l - poo) + . (20) 
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Obviously, this combination is invariant under rotation around the y' axis. In fact, it is proportional to one of three 
eigenvalues of the matrix p (see, e.g., @) so that 



H 

Poo 



1 H 

2 Pi-i 



K 

Poo 



1 K 

2 Pi-i . 



(21) 



where p^m' an< ^ Pmm' are the correlation coefficients calculated in the H- and i^-system. 

Equation (| 19() is the basic equation that will be used in the rest of the paper to evaluate the angular distributions. 
Its structure is independent of the particular frame chosen for the description of the A decay, since it is fixed by the 
A spin and the parity conservation condition j6]). The coefficients p m m'{&) are determined by the spin-parity of the 
resonance R. Therefore, analysis of the experimental angular dependence should enable one to get information about 
the resonance quantum numbers. 

A good test of the production mechanisms might be the distribution Fy(0; 9, (j)) at = 0, where the A decay is 
observed at forward direction in coincidence with r\ moving in the opposite direction along the beam axis. Then as 
follows from Eq. (fl"5|) the matrix p becomes diagonal and the cj> dependence in W(Q;9, (f>) disappears (at O = the 
cross section is obviously invariant under rotation around the Z axis) 



W(Q = 0;9,4>) = 

= ^ (am)(0) cos 2 9 + ~(1 - poo(0)) sin 2 6>) 

= ^ (poo(0)(1 + 3 cos 2 9) + (2 - 3,000(0)) sin 2 0) 



(22) 



with poo(0) = Poo(0 = 0). Using m = ml = in (|15p we obtain the following formula for poo(0) 

2 



Poo(0) 



3(1 + ac) 



(23) 



with 



Cr 



_ 3 2(i+l-j); 



C 



J - 



L0i 



^3/2 



(24) 



One can readily see from Eqs. (|22ll to (|24|) that the first term in Eq. (|22|) . proportional to 1 + 3 cos 2 9, is provided 
exclusively by the A helicity p,A = 1/2 whereas the second term with sin 2 9 is due to p\ = 3/2. 

The expressions for f(9) — 2ttW(Q = 0;9,4>) are summarized in Table UTl The situation is particularly simple 
for J = 1/2. In this case f(9) does not depend on the electromagnetic part and its form is totally fixed by the A 
spin. It is identical to the angular distribution of pions through the A decay in single pion photoproduction. For 
higher resonances having J > 3/2, also the substates with the helicity p,A = ±3/2 become populated. As a result, the 

element poo and consequently the shape of f(9) depends on the ratio a = • particular, it is convex 

upwards (downwards) in the whole region — 1 < cos 9 < 1 if ac < 1 (> 1). According to the formulas in Table UTl the 
angular distribution for J* = 3/2 + and 5/2~ is convex upwards for not too low values of a. In the other two cases 
J n = 3/2" and 5/2 + we should observe quite a slight angular dependence with a convex shape downwards for a < 1 
and a < 3/2, respectively. Taking, e.g., A 3/2 /A 1/2 = 0.81 for D 33 (1700) from the PDG 0], we obtain f{9) w 3 + cos 2 9. 

Besides relatively simple formalism (Eqs. (f2"2"j) to (fM]) ) the measurement at = has the advantage that at very 
forward angles the overlap between ttN and rjN states becomes minimal. It is especially important at low energies, 
where the restricted phase space does not allow the particles, e.g., r\ and A to escape the interaction region before the 
A decays. Therefore, this method is a possibility to naturally reduce the corrections appearing when the A decay is 
influenced by the presence of the 77 meson. 

As a next step we consider the distribution over the angles 9 and <p of the pion momentum in the nN rest frame. 
First, we introduce a new distribution function defined as 



W(fl) = ^-f \t mf x(Q;n)\ 2 sin0rf0, 



w(n)dn = 1 , 



with the normalization constant 



N = 



( A V2 



^3/2 



(25) 



(26) 
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It is clear that the general structure of W(Q; f2) (|19p also holds for VF(fi) so that we can immediately write 

in in 

mm' 

= (poo cos 2 9 + ^(1 — poo) sin 2 8 — \f21Zepio sin 28 cos 2(f) — sin 2 # cos 2^) . (27) 
An V 2 / 

Due to parity conservation the third term in the brackets proportional to sin 28 should vanish, because it changes 
sign under the transformation 8 — > n — 8, <f> — * (f> + n. Therefore pio = and the matrix p has only two independent 
elements po = poo and p\ = f>\-\ and is of the form 

|(1 -Po) px 
P = | po | • (28) 
Pi |(1 -po) 

It is instructive to consider the eigenvalues a, (3, and 7 of the matrix p 

«= Tj(l-Po) -Pi, /3 = i(l-p )+/>i, 7 = Po, 

a + + 7 = 1. (29) 

The eigenvalue /3 is proportional to the combination (|20|) integrated over 0. From (f29|) and positive definition of the 
matrix p the following restrictions hold in any reference system 

0< po <1, < |(1 -Po)- (30) 

These inequalities can be useful in reconstructing the matrix p (|28[) from the experimental angular distributions on 8 
and 4>. 

Using Eq. ((4]) we obtain the following formula for the correlation coefficients in the helicity frame 

n H = — WJ^' r J , n H = 1 f31) 

Hmm' 2]\[ / ; / i ni(A,mm' u >.w ' / , Pmm A ) V OJ -; 

m^A /i/i' 



where 



^mm' — n J i x ' r<2^ \ aR\ 2 f^oA 



2J+1 



2 



7T 

/ <(©)<'(©) sme<ie : 





(33) 



A,— — M A./^/i 

Substituting to = to' = and to = — m' = 1 into (|3 1 [) it is straightforward to obtain the values of po and pi for an 
individual partial wave i?(J 7r ). In particular, the parameter po is of very simple form 

3 \2J+l) V°io§i 



f$ = P& = o l^T-hr) KoV) . ( 34 ) 



so that the diagonal elements of the matrix p (f2"B")) are totally independent of the electromagnetic part and are fixed 
only by the spin-parity of the resonance R. For pi we have 



jLtjLt— 2 



X U L0| M iO| M -2 L/ llIm / °l-lim / ' V J0 J 
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or using symmetry properties of the Clebsh-Gordan coefficients and the relation (|33| for the coefficients c J 



- ± - 

Since c| 3_j_ = 0, the resonances with J 71 " = i do not exhibit (f> dependence. Furthermore, cf 3 _i = 0, so that for 

2 ' 2 2 2 ' 2 2 

J = 5/2 the second term in the brackets, corresponding to the initial helicity state A = 3/2, does not contribute. This 
leads to the fact that the sign of p\ for J 71 " = 5/2 ± does not depend on the parameter a and therefore the slope of the 
corresponding (^-distribution at cf> — ► (or <j) — ► 7r) is model independent. 



In the canonical system 



Pmm ! 



2iV 

m/A M L M' 



where 



„M L M' L _ n ]\ n J\ r ,\M<± n \M'^ \ n R aR\ 2 (oo\ 



or i i 

4m;, = — / 4f i0 (©) djff o(0) ^ 6 d6 . (39) 



o 



The elements Pq and are then given by 

_ 2A 7 " ^— ^ D m/A,0OJ 



m/A v 



Using the last equations one can easily find 

2 

3(1+ a) 



P? = ^-^ F [fCff* 3 )\a(c J r ! 3 ) 2 ), (41) 



m/=±i 



It follows from (j4"2"|) that the amplitude A?, 2 can contribute to the nondiagonal term p± only if L > 2. In our case 
this is only ,V = 5/2" (see TableHJ). 

The values p^ H and p^ H for all six states are given in Table HTT1 One can see that the combination po — 2/?i 
does not change if we turn from the K- to the iJ-system. This fact becomes trivial if we notice that the equality 
Pq — 2pi = p^ — 2pf" immediately follows from the Eq. (|2Tj) after integration over 9. 

Now let us consider the distribution over cos#. Its structure is similar to that in the previously discussed case (f2"2")l 

W{9) = J W(8, <jy) dtp 

= I (po(l + 3 cos 2 9) + (2 - 3p ) sin 2 6>) . (43) 

In analogy to W(Q — 0; 0) the first and the second terms in Eq. ([43]) . evaluated in the -ff-system, are related to the A 
helicities \pa\ — 1/2 and \pa\ = 3/2, respectively. It is important that in the 7i-system the element p$ is independent 
of a and is totally determined by the quantum numbers of R (see Table HIT)) . Furthermore, resonances with different 
angular momentum and parity contribute incoherently to p H . Therefore, the function W(9) calculated in the helicity 
frame is especially effective as a tool to identify R. 
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TABLE III: Coefficients po = Poo and pi = pi-i in Eg. (127fl (see also ((43} and (144 [l l calculated in the H- and the if-system. 



J' r (Z/2T2j) 



f (S31) | + (f3l) 



(-Da 



1 + 
2 



(£3 



15 



// 
Pi 



1 1—0 

"6 l+o 



1 1 
7 l+o 



1 1 
' 5 l+o 



2 1 

3 1+a 



_2_ 3+2a 
15 1+a 



1 3+2a 
5 1 + a 



pf 



_4 1_ 

15 1 + a 



70 1 + a 



J 1_ 

10 1 + a 



As for the (j> dependence of the cross section, its structure follows from the general expression (f27j) 

W{4>)= [w(6,<p)sm6d6 = —(l-2p 1 cos2<f>), (44) 
J ^ 

where p\ in H- and AT-systems are given by the Eqs. and (142[) . As already pointed out, the <f> dependence for ^ 

is trivial in the helicity system. Furthermore, for the highest resonances with ,V — | , although the amplitude of the 

oscillations of W depends on the parameter a the character of its convexity in the region < <p < ir is independent 
of a. In the -ftT-system, the 3/2-helicity amplitude contributes to pi only for J T = 5/2 _ . In the whole, as one can see 
from Table IIII1 the sign of p\ is fixed only by the spin-parity of the resonance and the character of the (j> distribution 
is model independent. 



III. ISOBAR MODEL FOR ~/N -> n°r]N. DISCUSSION OF THE RESULTS. 

The expressions presented in the previous section relate to the ideal situation where the amplitude for 7^ — > ir°T]N 
is dominated by the single diagram, Fig. [ljg). The natural question arises: what is the influence of the 7r^V* channel? 
Inserting the corresponding diagram, Fig. [ljh), into the formulas above will certainly make them much lengthier 
and less symmetric. Therefore we consider the problem numerically and discuss in this section the influence of the 
ttN* production on the results predicted by the formalism of Sect.[nJ For further study we need a model describing 
photoexcitation of the state i?(J T ) and its decay according to the schemes (a) and (b) in Eq. ©. Here we adopt a 
typical isobar model along the line used for double pion production (see, e.g., p, 0, §)■ 

The amplitude used in the calculation is a sum of the eight terms corresponding to the diagrams in Fig. [1] 

J 1 = J"(a-/) _|_ J^ia) _|_ y(>i)_ (45) 

The first six graphs (a-f) form the background. We neglect the diagrams with ijNN coupling due to its weakness. 
Therefore, the main model parameters are the partial widths of N* . We use 45 % for both r]N and ir N modes 
and 10% for that of 7r7rA^. The total width is equal to Tn* = 150 MeV. As already mentioned, the background 
mechanisms provide only a small fraction of the observed cross section for 773 — * r)ir°p, and this fact is considered 
as a key indication that the reaction mainly proceeds through resonance excitation. The corresponding amplitudes, 
depicted in Fig. HKg-h), read 

T(»> = T lN ^ A G A F A ^ N , (46) 



with G A and Gat* standing for the A and Af* propagators. 



(47) 
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For each resonance i?(J 7r ) we used a simple nonrelativistic Breit-Wigner ansatz with an energy-dependent width 

r«(w) 

{X\T jN ^ vA \q, M A ) = A* G R (JA I^r^aI?, M a ) , 
(AIT^-^jv. |<f, Mat.) = A$ G R (JA jF^,^, |g, M N ,) , 

G R = (w - M R + ~r R (W)\ , A = 1/2, 3/2. (48) 
The vertices Fr_> x (x 6 {i]A, irN*}) in Eq. (|48[) were taken in the phenomenological form 

(JM\Fr^ a \ q, M A ) - C££ £ 3 Ma F l * Ml (q), L = L( V A) , (49) 

(JM|F R ^. I £ Mv*) = f RwN , £^ , Mjv , r L V (9)> i = M**r) . (5°) 

where L(r]A) and L(irN*) are given in Table [H As already noted, in (|49l we assume only the lower of two possible 
values of the angular momentum L(j]A). In both channels (|4"9"|) and the finite width of the A and N* isobars 
was taken into account, what is important for the low energies considered here. Significant contribution to the width 
of each resonance is assumed to come from the tvttN mode (see Eq. (|54[1 ). The corresponding energy dependence was 
taken in a simple form 

T wwN (W) ~(W- M N ~ 2m^) Q(W - M N - 2m 7r ). (51) 

Since our calculation relates to the region of low kinetic energies we use the nonrelativistic formalism for A and N* 
states. Therefore, we do not touch upon such a complication as off-shell ambiguity in the Fa^ttN vertex in Eq. (|46j) 
appearing in the relativistic treatment of the spin 3/2 field (see, e.g., Ref. Q). 

Taking the R — > irN* transition in the phenomenological form (|50[) we have essentially simplified the problem in 
comparison to Ref. where the D33 — > nN* decay is calculated microscopically. Within the approach of Ref. 
the process 7iV — > U33 — > 77 A is treated as a driving mechanism producing the Tr°r]N state already at tree level. The 
coupling constant Jr v a entering the vertex R — > r/A is taken from the analysis of Ref. (Io| . Then the production of 
the ttN* state proceeds as a series of interactions D33 — * r/A + — > n a rjp — ► tt°N* , taken up to the first order in the 
corresponding two-body scattering matrices. In our case, the constant /r^n* hi (f50|) is real and its absolute value is 
fixed by the unitarity condition which for the Breit-Wigner resonance reads 

T^AT* = Tr — T n N — T„„ff — T, ; A ■ (52) 

The model should be reasonably good at least close to the resonance position W « Mr. With distance from this 
point the energy dependence dictated by the Breit-Wigner ansatz (|4"5|) could differ from the one obtained within 
the microscopic approach. Clearly, the most unambiguous treatment would be a three-body calculation including all 
coupled channels rnrN, ttN and irrjN. 

Close to the ir-q production threshold, the functions T vA (W) and T^n* (W) are mainly determined by the centrifugal 
barrier effect resulting in T x ~ q% +1 , where x € {77A, irN*} and the orbital momenta L(x) are collected in TableQ] 
We can expect that in the low energy region the relative fraction of the ttN* channel is important in the ^ channels 

since their decaying into 77 A requires higher values of L(r]A). As already noted, the | state might be a well candidate 
to explain the rapid rise of the cross section in the threshold region. Another state producing s-waves in the irrjN 
system is Any appreciable amount of other states is less likely since their decay into nriN at low energies is 
suppressed by the centrifugal barrier. 

In Fig.[3] we show an example of the total cross section calculated with J* = | . The calculation demonstrates 
a strong dominance of the resonance mechanism over the background terms (dashed curve). This result agrees with 
that of Ref. @. 

For the helicity amplitudes of -D33 (1700) we used average values [|J 

A 1/2 =0.104 GeV~ 1/2 , A 3/2 = 0.085 GeV~ 1/2 , (53) 



and for the mass and widths 



Mr = 1.72 GeV, rR = 300MeV, -^=20%, 

1 R 

%^ = 7 3%, ^ = 5%, -^-=2%, (54) 

I R 1 R lii 
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FIG. 3: Total cross section for jp — > ir°rip. The dashed curve is the background contribution (diagrams (a) to (f) in Fig.[T|). 



where all T's are taken at W = Mr. The last three values were chosen simply by adjusting the resulting total cross 
section for — > n°r)p to the data of Ref. [l[ . It is remarkable that fitting the cross section may require quite a small 
strength of £>33(1700) decay into the r/A and irN* channels. In order to get a feeling of possible variation of T v a and 
rVjv* one can use the formula for the total cross section at the resonance position 

a(M R ) = C T (2J+l) ^f 2 ^ N , (55) 

where all energy dependent quantities are calculated at W = M R . The coefficient Ct takes into account the isospin 
structure and T yN is the radiation decay width. Taking for £> 33 (1700): C T =4/9, M R = 1.7 GeV, T R = 300 MeV, 
and T jN /T R = 0.19% from '4\ we obtain 

« 8.5 • 10-V(M i? ) . (56) 

1 R 

Since around the total energy W = 1.7 GeV the cross section has a strong rise the ratio is very sensitive to Mr. 
For instance, for M R — 1.72 GeV we will have a(M R ) w 0.8 rib and T^^n — 7%T R . The maximum value a ~ 4 fib 1] 
gives about 34% for the total nrjN width. Then the inequality 

r„ A + iw < r^jv = 0.34 r R (57) 

holding for the constructive interference between rjA and TtN* configurations gives an upper limit for the sum of the 
partial decay widths in 77A and TtN*. 

In the following we discuss a general case in which a resonance R(J V ) produces mriN according to the two schemes 
in Eq. ([2]) and the background is totally neglected. For each resonance R{ J 7 *) we use the same parameters 

M R = 1.8 GeV , T R = 300 MeV , = 20% , 

Tn7rN - 60% , ^ = 20% , ^ = 10% . (58) 



r_R, T R r R 



For the ratio a (flOt we take 



» S (^) J =0.67. m 

Firstly, we show in Fig.|¥]the tcN invariant mass spectrum where the contributions of the final state configurations rjA 
and ttN* (schemes (a) and (b) in Eq. ^) are separately presented. It is interesting that the overlap of these states 
essentially differs in different partial waves. It is quite large in J v = § and I and less essential in other waves. 
Clearly, the character of the interference depends on the particular values of the orbital momenta L(r)A) and L(ttN*) 
as well as on the relative sign of the Rr/A and RirN* coupling constants. 
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FIG. 4: Distribution of the 7riV invariant mass calculated at the total jN cm. energy W = 1.8 GeV (corresponds to the photon 
lab energy about 1.255 GeV). The dashed and the dotted lines show the contribution from the r/A and nN* production. The 
solid line is the coherent sum of both channels. The results are presented in arbitrary units. 
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FIG. 5: Angular distributions f(8) — W(Q — 0; 8, (f>) of n mesons in the nN cm. system when the angle O is fixed to = 0. 
The dashed curve contains only the contribution of the driving 7/ A term in R — > TTrjN decay (the channel (a) in Eq. ((2J). In the 
dotted curve also the nN* channel ((6) in Eq. ([2])) is taken into account. The solid curve represents the symmetrized function 
fs{8) (see Eq. (|60[0 where in addition the invariant nN energy is restricted to the region M^n < 1-13 GeV. The total 7./V cm. 
energy is W = 1.8 GeV. 
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FIG. 6: Distribution W(9) (Eq. ((43])) over the polar pion angle in the nN rest frame calculated in the helicity system at W — 1.8 
GeV. Notation of the curves as in Fig. [5] 
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FIG. 7: Distribution nW(<f>) (Eq. (I44|l ') over the azimuthal pion angle in the ttN rest frame calculated in the helicity system at 
W — 1.8 GeV. Notation of the curves as in Fig. [5] 




FIG. 8: 7r°?7 invariant mass spectrum at W = 1.8 GeV for different transitions listed in Table U in arbitrary units. The results 
are obtained without taking the nN* channel into account. 



Of special importance for us is the contribution of the ttN* state in the region of the low invariant masses M v n. 
As we can see in Fig. @] it is important in i + and | states and can be neglected in other waves if sufficiently low 
values of M^n are considered. Therefore, to isolate the dynamics related to the scheme (a) in Eq. @ we only need 
to exclude the region with M^n larger than a certain value M® N depending on the overall reaction energy W . Then 
it is reasonable to assume that the A life time is sufficient to escape interaction with the r\ meson. In the waves with 
J OT = ^ and | the situation is more complicated at least at the energy W = 1.8 GeV considered here. Clearly the 
overlap between 77A and irN* should decrease with increasing W. 

In the series of figures [5] to [7] we present examples of angular distributions for ir°r] photoproduction by unpolarized 
photons on an unpolarized nucleon. The calculation is performed at W = 1.8 GeV. In all figures the dashed curve 
represent the distribution, in which only the channel (a) in Eq. @ is taken into account. Their form is described by 
the analytic expressions obtained in Sect.|TTl Addition of the scheme (b) gives the dotted curve and the solid curve 
is obtained after cutting off the kinematical region with M t jv > M®n = 1.13 GeV. As expected, the interference 
between the channels (a) and (b) tends to distort simple angular dependence, obtained if only the scheme (a) is used. 
After eliminating the energy region in which both mechanisms strongly overlap, we bring the calculation back to 
qualitative agreement with the results shown by the dashed curves. This effect is observed in all waves except for 
J n = i + and I in accord with our notion about strong overlap of the rj A and nN* configurations in these states. In 
the whole, using the above procedure the qualitative features of the resulting angular distributions are in agreement 
with the simplified calculations, in which the ttN* channel is neglected. 
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The solid curve in Fig. [5] contains symmetrization with respect to cos 9 = 0, i.e. 

f(6)^f s (6) = ±(f(B)+f(ir-0j). (60) 

In general, after addition of the channel nN* the resulting angular dependence yields certain forward-backward 
asymmetry, which can make the analysis more complicated. This effect is removed after the artificial symmetrization 

The ^-distribution in Fig. [5] is quite similar to the one in Fig. [5] As already noted the convex up and down form of 
W(9) related to sin 2 9 and (1 + 3 cos 2 9) terms in Eq. |J52J is provided by the A helicities |/ua| = 3/2 and |^a| = 1/2. 
Flat distribution for J 71 = | is due to the fact that in this instance (after integration over cos 0) we have an even 
mixture of both helicity states so that the sum of |(1 + 3 cos 2 9) and sin 2 9 gives a constant value. In other words, in 
the ii-system the longitudinal A polarization (averaged over the region — 1 < cos 9 < 1) is zero if J n = 3/2 - . In the 
case of a small admixture of the 7riV* background under the A peak we will have an interference term proportional 
to cos 9 caused by different parities of N* and A. This effect is clearly seen in Fig. [6] (dotted line in the panel for 
J v = | ). In other cases addition of the irN* channel leads to more complicated form of the function Wiff). 

Knowledge of the angular distribution W(9) in the helicity system is of key importance for understanding the 
spectrum of the ir-q pairs. Indeed, for each mass M^at the value of cos 9 (9 = 8*) is determined by M nrj [ll[ according 
to 



cos 8 



1 



AWAU N q*q v 



(M' N - M Z N + mi){W z - M' N - m*) - 2M^ N (M^ - m* - mj)) , (61) 



Therefore, the knowledge of W(cos8) at fixed M^jv immediately provides the Dalitz plot distribution 
d a/dM^^dM^Tj, and for the ttt) spectrum we have 

d(J = f d 2 a MnN , M (f , 2] 

dM% n J dM nN d cos 8 2Wq*q v * N ' 1 ' 

In other words, if there are no nearby resonances in the irr) system (as in our case), the structure of the Dalitz plot 
(M^jv, M^ri) is totally determined by the quantum numbers of the resonance R related to its decay into the 77 A and 
the ttN* channels. In Fig.[8]we present the spectrum da/dM 2 N given by different states i?(J 7r ). As we can see, apart 
from the boundary of the allowed kinematical region where da/dM 2 v —> the spectrum qualitatively reproduces 
the shape of the angular distribution in the helicity system (dashed curve in Fig. [5]). Thus the nrj mass distribution 
should be sensitive to the quantum numbers of the resonance R. Again for J 77 ' = l/2 ± the spectrum, having a visible 
minimum in the middle part is independent of the electromagnetic properties of the resonance. For other resonances 
it depends on the parameter a = (Ay 2 /Ay 2 ) 2 . 

As for the invariant mass distribution in other two-body subsystems, the corresponding measurements can hardly 
give useful information. As an example, we can take da/dM^N shown in Fig.|4j First of all, the general structure is 
quite insensitive to the choice of R(J^). If we change from one resonance to another, in the main only the position of 
the maximum is shifted. Obviously, this shift is explained by the barrier effects. Namely, since da/dM^^ ~ da/dq v , 

for low values of q n the spectrum is proportional to , where L(rjA) is the angular momentum of the decay 

R — > 77A (see Table IJ). With increasing L(rjA) the centrifugal barrier factor tends to suppress the cross section at 
low q v (large M^jv), resulting in shifting the maximum to higher values of q n (lower M^n)- This trivial effect is what 
we mainly observe in Fig.|4j Furthermore, at higher energies the shape of the spectrum around M^n = Ma will be 
governed by the nN energy distribution in the A region, so that the values of da/dM^ N is mainly determined by the 
form of the A peak. In this connection, investigation of da/dM^N is not of any use to get additional information on 
the reaction mechanism. 



IV. CONCLUSION 



We have discussed some details of a phenomenological analysis of ^N — > ir°r]N aimed at identifying the dominant 
mechanisms of this reaction. This analysis assumes that the part of the amplitude corresponding to a given resonance 
R is sufficiently large, so that other partial waves and the background can be neglected. This assumption seems to 
be justified by direct calculation of the most important Born diagrams (see Fig.[T]). Our results show that even in the 
absence of the polarization data some interesting properties of the reaction can be found. Our main focus is on the 
angular dependence as a test of the mechanism responsible for ir r] photoproduction. Here we summarize the most 
important qualitative features of different types of the angular distributions considered in the main part of the paper. 
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1) The distribution W(Q) over the polar angle of the irN system in the jN cm. frame (TableHH). In the simplest 

case of J Tr — 5 we have W(Q) = const. Generally, the distribution on depends on the parameter a = (A 3 / 2 / 'Ai/2) 2 \ 
which hampers the analysis if the electromagnetic properties of the resonance are poorly known. Indirect clues to the 
spin-parity of the resonance can be gained from the complexity of W(Q). For a state with 77A angular momentum L 
the complexity of the VF(0) is fixed to 2L for J > |. 

2) The distribution over cos 9 at fixed = (Fig.[5|). If we select 77 mesons produced opposite in the direction to 
the photon beam, a simple formula (|23[) can be obtained. In the general case the reaction mechanism is such that 
the intermediate A states with both helicities \pa\ = 1/2 and 3/2 are populated. As a result, the decay angular 
distribution of A differs from the simple form 1 + 3 cos 2 9, peculiar for single pion photoproduction. The weight of 
each /ia configuration depends on R( J 71 ') so that the method should allow the quantum numbers of R to be extracted 
from the measurements. As in the W(Q) case, the particular form of the distribution in the states with J* = | and 
I depends on the value of the parameter a, Eq. (jTUJ) . However, in a wide range of a, a > | for ,V = | + and a > i 
for J 7r = | , the corresponding distributions are not very sensitive to a. 

3) The distributions W(9) and W((f>) over the polar and azimuthal pion angles in the ttN rest frame. For convenience, 
we express the functions W(9) and W((f>) in terms of the A decay correlation coefficients p mm >. Their values can be 
determined by fitting the analytic expressions for W to the experimental data. 

For W(9) the helicity system (Fig. [6]) seems to be more useful for the partial wave analysis since the decay corre- 
lation coefficients p!^ m , and the corresponding angular distributions are independent of a. Furthermore, there is no 
interference between partial waves with different spin-parity J'*, so that resonances contribute incoherently. This is 
especially important in the situation of strongly overlapping states, which is quite typical in the second and the third 
resonance region. However, there is a sensitivity of the angular distributions to even small admixtures of the itN* 
channel. Furthermore, in this case there are qualitative difficulties in distinguishing between the states belonging to 
the two groups J* = l/2 ± ,5/2+ and J n = 3/2^5/2". 

The distribution W((f>) in the iJ-systcm for J 71 " = | seems to be sensitive to the parameter a, Eq. (jTUJ) . This 
shortcoming is however partially avoided in the canonical frame where the sign of the ^-dependent term does not 
depend on a. 

It has also been shown that if the correlation coefficients p~Q are fitted to the data, the spectrum da/dM^N 
does not provide additional information, since the Dalitz plot is immediately obtained if the distribution over cos 9 in 
the ii-system is known. This result is a trivial consequence of a linear relation between the cosine of the pion decay 
angle and M% n (see Eq. (J6T|) ) . 

In general, our results demonstrate that different assumptions about the spin-parity of the dominating partial wave 
lead to different predictions regarding angular distributions, so that each state R shows its own signature. Model 
independence of some of these signals as predicted by the present analysis is the major motivation for proposing our 
method. Additional important information can be obtained from the polarization experiments. 

Finally, we would like to note, that the presented calculations are related to the case when most of the Tr°rjN 
configurations are produced through the R — > rjA decay. On the other hand, rather low percentage of r/A and nN* 
channels in Eq. (f54"|) indicates that the the measured cross section [l[ may be accounted for with quite a small fraction 
of 77A and 7riV* channels in the total resonance width. This observation might be a reason to doubt the necessity 
of introducing direct RrjA coupling to explain ir°r] production rate. The needed strength can be provided by the 
R — > 7rA and R —> pN decays followed by irN — * r/N rescattering. In this case we do not need to restrict ourself to 
the resonance states with the isospin T — 3/2. Clearly, another assumption about the reaction mechanism, as to the 
one in which the ttN* dominates, will lead to angular distributions different from those presented here. 
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